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A decimation method is applied to the tight binding model 
describing the two dimensional electron gas with next nearest 
neighbor interaction in the presence of an inverse golden mean 
magnetic flux. The critical phase with fractal spectrum and 
wave function exists in a finite window in two-dimensional 
parameter space. There are special points on the boundaries 
as well as inside the critical phase where the renormalization 
flow exhibits higher order limit cycles. Our numerical results 
suggest that most of the critical phase is characterized by a 
strange attractor of the renormalization equations. 
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I. INTRODUCTION 

The two-dimensional electron gas with irrational mag- 
netic flux is a well-known paradigm in the study of sys- 
tems with two competing periodicities. The magnetic 
field results in reducing the problem to a one-dimensional 
tight binding model (TBM) |^ known as the Harper 
equation. Q The Harper equation exhibits both extended 
(E) and localized (L) states. At the onset of transi- 
tion corresponding to a periodic potential with square 
symmetry, the states are critical (C) with fractal spectra 
and wave functions. The scaling properties of the devil 
staircase spectra and the wave functions have been stud- 
ied extensively using various renormalization group (RG) 
methods. 1, §, § 

Recently, it was pointed out that the Harper equation 
also describes the isotropic XY quantum spin model in 
a modulating magnetic field of periodicity incommensu- 
rate with the periodicity of the lattice. ||] It was shown 
that the presence of anisotropy in spin space fattened 
the critical point of the Harper equation resulting in a 
phase diagram where E, L, and C phases all existed in 
a finite measure parameter interval. The existence of a 
fat C phase provided a new scenario for the breakdown 
of analyticity in incommensurate systems. Furthermore, 
based on numerical results obtained using a new decima- 
tion method 0, ^ it was argued that the fat C phase 
was described by four distinct universality classes char- 
acterized by limit cycles of the RG flow. 

The fat C phase was also reported recently in the TBM 
describing Bloch electrons moving on a tight binding 



square lattice where the coupling to next nearest neigh- 
bor (NNN) sites exceeded a certain threshold value com- 
pared to the nearest neighbor (NN) coupling. Using 
various analytical and numerical tools to study the scal- 
ing properties of the fractal eigenspectrum, the universal- 
ity classes of the C phase were investigated. The E-C and 
C-L transition lines were conclusively shown to define a 
new universality class, different from the E-L transition 
line (which belonged to the Harper class). However, the 
analysis in the interior of the C phase was rather incon- 
clusive. 

In this paper, we apply decimation methods to the 
TBM describing Bloch electrons in a field with a NNN 
interaction. Unlike the previous study ]lO| which inves- 
tigated the scaling properties of the fractal eigenvalues, 
we study the scaling properties of the fractal eigenstates. 
Our study, in agreement with the previous work ]l"o|, 
shows that the E-C transition line defines a new univer- 
sality class while the E-L transition line belongs to the 
universality class of the Harper equation. In addition, we 
show that there exist special points in the interior as well 
as on the boundary of the C phase which correspond to 
many new universality classes. They are associated with 
higher order limit cycles of the RG equations. Our de- 
tailed analysis suggests that the rest of the interior of 
the C phase is described by a unique universal strange 
attractor of the RG flow. 

In Section II, we describe the model and briefly review 
previously known features of the phase diagram of the 
TBM describing Bloch electrons with NNN interaction 
and make a comparision with the quantum spin prob- 
lem. The decimation scheme is reviewed in Section HI. 
Section IV concentrates on higher order symmetric limit 
cycles, and in Section V we describe the phenomenon of 
shifted symmetry and doubling of the period of a limit 
cycle. Evidence for a strange set is provided in Section 
VI. Finally, Section VII contains a summary of our con- 
clusions. 



II. TIGHT BINDING MODEL FOR BLOCH 
ELECTRONS IN A FIELD 

The TBM describing Bloch electrons on a square lat- 
tice in a magnetic field with both NN hoppings ta and tb 
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and NNN hoppings tat = <ah is 0] 

{ta + 2tabCOs[2TT{a{i +^)+ <l>)]}^^+l 

+{ta + 2tabCos[2n{aii ~ ^) + ?!')]}V'i-i 

+2tbCos[2TT{ai + Eip^ (1) 

Here, a is the magnetic flux which we choose to be the 
inverse golden mean a = {Vb — l)/2. This TBM was 
studied in detail in ref. . Fig. 1 shows the phase dia- 
gram of the model in the space of the parameters A = |^ 
and a — 2y^ . The Harper equation corresponds to the 
limit a = where the NNN coupling term is zero. The 
phase diagram was obtained using analytical meth- 
ods to obtain the scaling behavior of the total band width 
(TBW) and the Lyapunov exponents and carrying out a 
numerical multifractal analysis. The lines AC (E-C tran- 
sition) and CE (C-L transition) were found to be bicrit- 
ical, i.e. the TBW scaled with the system size with the 
exponent 5 — 2. This was in contrast with the critical line 
BC separating the E and L phases where the exponent 
was known to be unity. Within the region bounded by 
the lines AC and CE and the a-axis, where the NNN cou- 
pling dominated, the multifractal analysis did not lead 
to conclusive results on the universality. This was due to 
lack in convergence of the f{a) curve with the size of the 
system. Furthermore, their numerical calculation of the 
TBW was complicated by oscillatory terms superimposed 
on the power law. However, this regime was conjectured 
to be critical. 

The existence of the fat C phase implies that the pres- 
ence of the NNN coupling in the Bloch electron problem 
introduces new universality classes. The same does not 
happen if the cosine term in the Harper equation is just 
replaced by more generic periodic functions. 

Eq. (1) involving both diagonal and off-diagonal disor- 
der bears some resemblance to the TBM describing quasi- 
particle fermion excitations in a quantum XY spin chain. 
1^ Unlike Harper, both the anisotropic XY spin chain 
and Eq. (1) exhibit C phase in a finite parameter inter- 
val. However, in the TBM (1) the fat C phase is observed 
beyond a critical value of the NNN hopping whereas in 
the spin model the fat C phase can be seen even with 
infinitesimal spin space anisotropy. 

Motivated by the success of our decimation scheme to 
confirm the existence of a fat C phase and obtaining the 
universal behavior in the quantum spin chain Q , we now 
apply the method to Eq. (1). Our main focuss is to 
obtain the universality classes of the fat C phase. 

III. DECIMATION SCHEME 

Our decimation approach describes the scaling prop- 
erties of the wave functions for a specific value of energy. 



Although all quantum states are fractal , it is more use- 
ful to study states at the band edges. This is because 
the self-similar behavior is usually observed only for the 
minimum and maximum energy states and also for the 
band center if i? = is an eigenenergy. In order to get 
the overall picture, it is sufficient to consider only one 
of these states and we will focuss on the quantum state 
corresponding to Emin ■ 

In RG analysis, in addition to fixing the quantum state 
, one has to also fix the phase factor (j> in Eq. (1). It has 
been pointed out in the previous studies |^ , , the wave 
function tpi obtained by iterating the TBM diverges un- 
less the phase factor cf) is tuned to some critical value. 
In the Harper model, the critical value of the phase fac- 
tor is ^ for the negative band edge. For this value, the 
main peak is centrally located and the wave function is 
symmetric about i = 0. In the study of the quantum 
spin model 0, the phase factor had to be varied contin- 
uously in the fat C phase so that the main peak could be 
centrally located and the resulting wave function became 
bounded. Determination of this critical phase factor was 
essential in order to find the RG limit cycles and to com- 
pute the universal scaling ratios. In general the phase 
factor (f> for obtaining symmetric wave function need not 
be identical to the phase factor resulting in bounded wave 
functions. 

We consider an infinite lattice which extends in both 
positive and negative directions from the i — site. In 
the decimation scheme, all sites except those labelled by 
positive as well as negative Fibonacci numbers are deci- 
mated. The resulting TBM connecting the wave function 
tp at two neighboring Fibonacci sites can be written as 

^{i + F„+i) = c+(z)V(i + F„) + d+«V(i) (2) 
^Pii - Fn+i) = c-{iWt - Fn) + d-iiWi). (3) 

The index n above refers to the level of decimation. 

For the Harper model it suffices to define only one set 
of the "decimation functions" c„(i) and dn{i) because of 
the symmetry of the wave function about i = 0. This 
implies that the scaling ratios are the same on the pos- 
itive and negative side. This was also asymptotically 
true in the quantum spin model everywhere else except 
along the C-L transition line 0|, where the wave func- 
tions were asymmetric about i = resulting in vanishing 
scaling ratios on one side, fl^ The reason why we now 
have to introduce separate decimation functions for the 
positive (+) and negative (— ) side is that sometimes the 
asymptotic (n — oo) c+(0),d+(0) appear to be shifted 
compared to c~(0), d~(0). This type of shifted symme- 
try turns out to be helpful in locating higher order limit 
cycles and is discussed in detail in Section V. 

Using the defining property of the Fibonacci numbers, 
Fn+i = F„ + Fn-1, the following recursion relations are 
obtained for c„ and dn (we will omit the 4-, — indices if 
the equations do not depend upon them) 0: 
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c„+i(i) = c„(z + Fn)Cn-l{i + Fn) - d~'^ {i)dn+iii) (4) 
dn+i{i) = -dn{i)[dn{i + Fn) + 

Cn{i + Fn)dn-l{i + Fn)]c-\i). (5) 

For a fixed i, the above coupled equations for the deci- 
mation functions define a RG fLow which asymptotically 
(n — > oo) converge on an attractor. In our earlier studies, 
the E, C, and L phases were distinguished by the distinc- 
tions in the attractors of the RG flow. In the Harper as 
well as in the quantum spin case, the C phase was char- 
acterized by a nontrivial asymptotic p-cycle at the band 
edges with p equal to 3 or 6. 

The existence of a nontrivial p-cycle for the decimation 
functions often implies that the wave function is neither 
extended nor localized and exhibits the self-similarity de- 
scribed by 

V'(^)«^(K^ + l/2]) (6) 

where [ ] denotes the integer part. The p-cycle of the 
self-similar bounded wave function can be used to define 
the universal scaling ratios 

= lim |V(Fp„+,)/V(0)|; J = 0, ...,p- 1. (7) 

n— *oo 

This equation describes the decay of the wave func- 
tion with respect to the central peak. A well-defined 
limit Q exists for an integer p for which asymptotically 
ip{Fn+p) ~ ip{Fn). For an even p, it often happens that 
|V'(-F'n+p/2)l ~ iV'l-Pri)! so that the above equation de- 
fines actually only p/2 different scaling ratios. Whenever 
any scaling ratio Qj takes a finite value between zero and 
unity, the wave function over the infinite system is nei- 
ther localized nor extended. In order to fully characterize 
the self-similarity of a wave function, an infinite number 
of scaling ratios have to be defined. 0] However, using 
the ones defined above one can already separate different 
universality classes from each other. 

The numerics to demonstrate a p-cycle was rather chal- 
lenging for the case where E = Q was not an eigenenergy. 
This was because at many points in the phase diagram, 
the energy was required with 16 digits precision (ma- 
chine double precision) in order to see the asymptotic 
cycle with two or more digits of precision. Even for tridi- 
agonal matrices, we were able to determine the energy 
only up to 12 digits. With this precision, the decimation 
equations could be iterated only about 16 times. The 
conjecture for the existence of a limit cycle provided a 
very efficient Newton method where the energy and the 
limit cycle were determined self-consistently. Diagonal- 
ization routines provided a good starting value of the 
energy which in principal could be improved to an ar- 
bitrary precision. At many points where the transients 
were rather long, the quadruple precision was used to 
confirm the existence of a limit cycle. 



IV. HIGHER ORDER SYMMETRIC LIMIT 
CYCLES 

Fig. 1 shows the E, C, and L phases of the model in the 
A — a space. The iteration of decimation equations shows 
that the EC critical line defining the boundary between 
the E and L phases is described by a 3-cycle of the RG 
fiow. On this line (with the exception of the point C) 
the decimation functions flow to the same limit cycle as 
for the critical Harper equation (point B). However, the 
point C is described by a different 3-cycle defining a new 
universality class. Table I compares the universal scaling 
ratios in these two universality classes. Fig. 2(a-b) show 
the wave function in these two cases. In both cases, a 
bounded wave function was obtained for (j) — 1/2 and 
the wave function as well as decimation functions were 
symmetric about the center of the lattice: i.e. c"*" = c~ 
and (i+ — d^ . 

The bicritical line AC was found to exhibit two new 
universality classes: At the point A, the decimation func- 
tions were found to asymptotically converge to a period- 
6 limit cycle for (p = 1/2. The regime bounded by the 
points A and C on the line AC (excluding the points A 
and C) appears to be described by a unique symmetric 
limit cycle of the period 12 for = 1/2. The 12-cycle is 
particularly clear for the middle point M (A = .5, a = 1). 
However, for other points on this line, the RG flow may 
exhibit long transients before settling on the limit cycle 
of the point M. This is illustrated in Fig. 3 for the point 
MM (A = .25, a = 1). The RG iterations begin close to 
the 6-cycle of the point A and then follow the 3-cycle of 
the point C for a while before approaching the 12-cycle 
of the point M . There is additional complexity involved 
because in order to observe the approach we had to shift 
the data for MM by 6 decimation levels. We used the 
quadruple precision to confirm that the MM limit cycle 
asymptotically converged to the limit cycle of M. 

Symmetric period-12 cycles were also observed at the 
points F, G, and H which fall on the line ta — tab- How- 
ever, the limit cycle at each point was different implying 
different scaling properties and hence different universal- 
ity classes. 

It should be noted that with the exception of the points 
B, C, and H, all the hmit cycles observed for 0=1/2 
resulted in asymptotically diverging wave functions (see 
Fig. 2). The decimation functions corresponding to di- 
verging wave functions cannot be used to obtain universal 
scaling ratios for the new universality classes. In section 
V, we determinine the bounded wave functions by tuning 
the phase factor. 

We also studied systematically other points on the lines 
CE and FH and points inside the C phase, but we did 
not see any evidence of limit cycles for (j) = 1/2- Since 
the numerics should be sufficient to show at least cycles 
of length 12 (unless there are very very long transients), 
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we can conjecture that there are no cycles of the order 12 
or shorter on these hnes with the phase 1/2. In section 
VI, we present some evidence that the rest of the C phase 
is described by a strange attractor of the RG flow. 

V. SHIFTED SYMMETRY AND PERIOD 
DOUBLING 

Fig. 4 shows the bounded wave functions obtained by 
tuning the phase factor to a certain critical value. Un- 
like the symmetric wave functions of section IV, these 
wave functions are not symmetric about i = 0. It turns 
out that this asymmetry is due to a constant shift be- 
tween the wave functions on the positive and negative 
sides. Numerical iteration of the TBM (1) shows that in 
addition to Eq. (7), the wave functions on the positive 
and negative sides are asymptotically related by 

^P{F„) « (8) 

Instead of the symmetric solution for the decimation 
functions corresponding to c+(0) — c^(0), d^{0) = d~{0) 
for all decimation levels n, we found the asymptotic equa- 
tions 

c+(0)«c-+,(0) (9) 
4(0) « d-+M, (10) 

i.e. asymptotically there was a shift of s levels between 
the positive and negative decimation functions. Morover, 
in above + and — could be interchanged with the same 
shift s, which imphcd that c+(0) « c~_|_g(0) ~ c^+2s(0)i 
i.e. the asymptotic period p = 2s. The shift s was found 
to be equal to the period of the symmetric limit cycle 
discussed in section IV. Therefore, the phenomenon of 
shifted symmetry resulted in doubling the period of a 
limit cycle for the decimation functions. 

The points F and G, which respectively fall on the 
intersection of the line ta = tab and the lines tb = and 
the self-dual line BD, exhibit the phenomenon of shifted 
symmetry with s = 12 (see Table II). This value implies 
the asymptotic cycle-length 24. 

The phenomenon of shifted symmetry is very crucial 
in locating the limit cycles of period 24. This is because 
with double precision arithcmatics, the RG equations can 
be iterated only about 24 times. Without having the 
shifted symmetry we could not have deduced that the 
asymptotic period is 24 as we could not go to big enough 
decimation levels to see the full cycle on one side only. 
This is shown explicitly in Table II. The same shift and 
period is observed for the critical phase at the point M 
in the middle of the line AC. At the point A the shift is 
6 and the asymptotic period therefore 12. 



VI. STRANGE SET OF THE 
RENORMALIZATION FLOW 

We explored the idea of describing the region bounded 
by the lines AC and CE by a strange set of the RG flow. 
Fig. 5 shows the two-dimensional projection of the at- 
tractor obtained by plotting an inverse decimation func- 
tion for two subsequent decimation levels in various parts 
of the phase diagram. Having Emin up to 12 digits, the 
RG equations are estimated to give correct decimation 
functions up to about 16 levels. Transients were taken 
into account by excluding the first six decimation levels 
from the data. It is interesting to note that the iteration 
of the decimation functions in three different parameter 
regimes namely the line FH, the region CHG (excluding 
the line CH), and the region ACGF appear to asymptot- 
ically converge on roughly the same set. Similar figures 
were obtained also on the line CG, GD and the lines EG 
and GH. The fact that different parts of the phase dia- 
gram are described by similar invariant sets makes us to 
exclude the possibility that the observed behavior is due 
to long transients. However, although the possibility of 
a very long limit cycle cannot be completely ruled out, 
we believe that the observed behavior suggests that the 
interior of the fat C phase ( excluding the special points 
which exhibit limit cycle ) is attracted to a unique invari- 
ant set of the RG equations. We conjecture that the set 
is a strange attractor. 

The iterates of the RG flow on the CE line seem to lie 
on the inner boundary of the invariant set correspond- 
ing to the interior of the C phase (see Fig. 5). In the 
previous studies |l^, the CE line defining the boundary 
of C and L phases was found to be bicritical. In anal- 
ogy with AC line, we would expect that the decimation 
functions on the CE line converge to a limit cycle of the 
order 24 for = 1/2. We did not see any evidence of this 
cycle. However, its existence can not be ruled out spe- 
cially in view of the possibility of long transients and the 
fact that even the symmetric limit cycle could be of order 
24. Therefore, the problem of determining the universal- 
ity class along the CE line describing the C-L transition 
remains open. 

VII. CONCLUSIONS 

In this paper, we demonstrate that our decimation 
scheme is an extremely useful tool to study general 
quasiperiodic TBMs. The Bloch electron with NNN in- 
teraction and the anisotropic quantum XY chain in a 
transverse field are two known examples where the C 
phase exists in a finite parameter range. In the spin 
problem, the C phase was characterized by four different 
limit cycles of the RG flow. |Q The present study shows 
that in the Bloch electron case the situation is lot more 
complex: in addition to six new universal limit cycles. 
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which correspond to self-similar wave functions, there is 
strong numerical evidence of a RG strange set. To best 
of our knowledge, this is the only example of a quasiperi- 
odic model where the golden mean incommensurability 
does not result in self-similar wave functions at the band 
edges. Even in the regime where no limit cycles exist, 
the RG scheme provides a clear distinction between the 
E, C, and L phases: in the E and L phases, the dec- 
imation functions are trivial, and in the C phase they 
assume finite non-trivial values. |Q It should be noted 
that in this regime where the C phase is not described 
by the limit cycles of the RG equations, the previous 
studies showed lack of convergence in /(a) curve. 

The general case of TBM (1) where the NNN inter- 
actions tab and are not equal provides an interesting 
limit of the Bloch electron on a triangular lattice. 
Our studies have shown that [0 the universality class 
of this model is related to the subconformal universality 
class of the Ising model. This establishes an interesting 
relationship between the anisotropic Bloch electron with 
NNN interaction and anisotropic quantum spin chains. 

In the quantum spin problem, the fattening of the 
C phase is due to the broken 0(2) symmetry in spin 
space which translates to the broken U(l) symmetry in 
the quasiparticle fermion Hamiltonian. For the Bloch 
electron, the existence of C phase is due to a NNN in- 
teraction. Although at this point we are unable to pin 
point the commonality between the symmetric breaking 
in these two problems, we believe that the origin of the 
fat C phase and new universality classes may be tied to 
certain broken symmetries of the models. 
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FIG. 2. (a-g) shows the wave function at the points B, C, 
A, M, F, G, and H corresponding to the phase factor <j) = 1/2 
such that the reflection symmetry about i = is preserved. 
This causes the wave function to diverge for the points A, M, 
F, and G while the wave function at the point B, C, and H 
is bounded. In these plots, the maximum value of the wave 
function is scaled to unity. 
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FIG. 3. (a) The inverse decimation function l/c„(0) vs. 
n {4> = 1/2) along the bicritical hne AC: The data for the 
point MM (A = .25, a — 1; shown by crosses), shifted by 6 
decimation levels, eventually follows the 12-cycle of the point 
M {X — .5,a = 1; solid line with small crosses showing the 
locations of the periodic orbit). 



TABLE II. The decimation functions for the point G 
{<f> = 1/3) at site i = showing the shifted symmetry with 
s = 12 and thus indirectly implying the limit cycle of length 
24. We see that c+(0) ^ c^±i2(0) and also d+(0) dn±i2(0). 



n c+(0) c-(0) d+(0) d-(0) 



FIG. 4. (a-d) show the absolute value of the wave func- 
tion at the points A = 1/4), M (0 = .3202185) , F 
{cf> = .2777...), and G (0 = 1/3) . The phase factor d> is 
chosen so that the main peak is centrally located resulting in 
a bounded wave function. 

FIG. 5. A two-dimensional projections of the inverse dec- 
imation functions inside the fat C phase. The data has been 
obtained by sampling on the line FH (a), in the region CHG 
(b), and (c) of the whole C phase which includes (a) and (b) 
and also the square ACGF and the points above it. The dark 
dots correspond to the data obtained along to the line CE 
(excluding the points C and H). 

TABLE I. The universal scaling ratios Q {j = 1, 3) at the 
point B (Harper) and at the point C (bicritical). 



j 0(5) Qic) 



6 


-0.817 


-2.477 


-3.499E-02 


0.109 


7 


2.852 


-6.304 


-1.257 


-3.952 


8 


-4.149 


5.509 


1.439 


-0.501 


9 


1.990 


1.176 


0.191 


-4.571E-02 


10 


-4.463 


-1.166 


-1.685 


0.302 


11 


2.733 


-4.122 


-1.168 


-1.808 


12 


-3.255 


42.235 


0.498 


-9.109 


13 


1.517 


1.559 


0.140 


-1.328E-02 


14 


4.982 


-6.744 


3.607 


0.250 


15 


-3.738 


-0.831 


-0.427 


-1.720E-02 


16 


-1.961 


-1.957 


0.467 


8.547E-02 


17 


-12.536 


9.281 


-6.087 


0.944 


18 


-2.507 


-0.809 


0.109 


-3.454E-2 


19 


-6.266 


2.874 


-3.942 


-1.26 


20 


5.532 


-4.132 


-0.502 


1.436 


21 


1.174 


1.996 


-4.567E-02 


0.191 


22 


-1.168 


-4.457 


0.303 


-1.684 


23 


-4.119 


2.738 


-1.807 


-1.168 


24 


42.288 


-3.256 


-9.106 


0.498 



(0, 2, 8, 34, 144, 610,...) 

1 (1, 3, 13, 55, 233,... ) 

2 (1, 5, 21, 89, 377,... ) 



0.2107 0.1712 

0.2107 0.2353 
0.2107 0.2387 
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